PN code chip time tracking with smart antenna

ABSTRACT

A pseudonoise code teaching loop (PNCTL) measures the difference in energy between a smart antenna array output that is despread using an early pseudonoise code and a late pseudonoise code and contains a chip time shifting current based on this difference. Using multiple elements of the array significantly improves performance.

BACKGROUND OF THE INVENTION Filed of the Invention

The present invention relates to spread spectrum communications generally, and more particularly to a method and apparatus for maintaining pseudo-noise code synchronization.

Smart antennas, consisting of multiple antenna elements and an adaptive algorithm, have received much attention recently because a smart antenna can suppress interfering signals if the direction of arrival angle (DOA) from the desired signal is different from that of interference, (Ayman F. Naguib, “Adaptive Antennas for CDMA Wireless Networks,” Ph.D. Dissertation Stanford University, August 1996 and Y. S. Song, H. M. Kwon, and B. J. Min, “Computationally Efficient Smart Antennas for CDMA Wireless Comm.”, IEEE Trans. on Vehicular Technology, Vol. 50, No. 6, November 2001.) Some encouraging field experiment results for smart antennas have been reported (M. Sawahashi, et al., “Enhanced Wireless Access Technologies and Experiments for W-CDMA Communications”, IEEE Personal Comm., pp. 6-16, December 2000). In a code division multiple access (CDMA) system, pseudonoise (PN) code synchronization is a difficult and critical process before data demodulation. The locally generated PN code should be synchronized with the received PN code for correct data demodulation. In general, PN code synchronization consists of two steps: PN code acquisition and PN code tracking. PN code acquisition is a coarse PN code synchronization to within a chip interval, and PN code tracking is a fine PN code synchronization (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995 and H. Meyr and G. Ascheid, “Synch. in Digital Communications”, John Wiley and Sons, 1990.) Most smart antenna studies have focused mainly on improving the bit error rate (BER) and capacity of a CDMA wireless communications system. Also, most of the time, perfect PN code synchronization has been assumed, or only one of the multiple antenna elements has been employed for PN code synchronization (Ayman F. Naguib, “Adaptive Antennas for CDMA Wireless Networks”, Ph.D. Dissertation Stanford University, August 1996; Y. S. Song, H. M. Kwon, and B. J. Min, “Computationally Efficient Smart Antennas for CDMA Wireless Comm.”, IEEE Trans. on Vehicular Technology, Vol. 50, No. 6, November 2001; Yoo S. Song, Hyuck M Kwon, and Byung J. Min, “Smart Antennas for 3G and Future Generation CDMA”, IEEE International Conference on Phased Array Systems and Technology, Dana Point, Calif., pp. 69-72, May 21-25, 2000) and M. Sawahashi, et al., “Enhanced Wireless Access Technologies and Experiments for W-CDMA Communications”, IEEE Personal Comm., pp. 6-16, December 2000.) This may be due to existing conventional antenna diversity schemes where antenna element spacing is large, e.g., ten times a wavelength, and arrived signals at different antenna elements are independently faded. This independent fading may make it difficult to coordinate the elements for the purpose of PN code synchronization. However, antenna element spacing in a smart antenna design for a CDMA system is usually assumed to be a half wavelength or less, and received signals at different antenna elements are the same except for phase differences due to the small spatial separation when the element response is identical. Calibration may be required, depending on the smart antenna algorithms, if the element response is not identical. If a normalized least mean square (NLMS) algorithm is used, then phase calibration is not necessary even though radio frequency (RF) mixers for each element have different phase responses (T. Chulajata, H. M. Kwon, and K. Y. Min, “Adaptive Antenna Array with No Phase Calibration for CDMA Reverse Link”, IEEE VTC'2000 Fall, Boston, Mass., Sep. 24-28, 2000). Also, mobiles may transmit signals with low power because of expected high smart antenna gains at a base station, and if only one element is employed for PN code synchronization, then the received signal-to-interference-plus-noise power ratio (SINR) at a base station may not be sufficient. As a result, a call drop may occur due to weak PN code tracking at a base station. It is more desirable for mobile users not to experience frequent call drops than to have high voice quality and fast PN code acquisition. These are among the motivations for a PN code chip time tracking loop (PNCTL) which can maintain its lock status for a longer time by utilizing all antenna elements instead of one element.

PN code acquisition using an antenna array has been addressed in other studies R. K. Madyastha and B. Aazhang, “Synch. and Detection of Spread Spectrum Signals in Multipath Channels Using Antenna Arrays”, IEEE MILCOM, pp. 1170-1174, 1995 and S. Kim and S. L. Miller, “Subspace-based Synch. Using Antenna Array for DS-CDMA Signals in Multipath Time Varying Fading Channels”, IEEE VTC Spring, pp. 547-551, Tokyo, May 2000) based on signal subspace techniques. However, many of these are difficult to implement because they require enormous computational loads, and also they assume short PN codes. The current and near future CDMA mobile applications e.g., the second generation (2G) and third (3G) CDMA systems (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995 and Third Generation Partnership Project; “Spreading and Modulation (FDD),” 3GPP Technical Specification, TS25.213, v3.2.0, March 2000) have employed long PN codes. Also, a practical PN code acquisition scheme employing all smart antenna elements has been studied in Bing Wang and Hyuck M. Kwon, “PN Code Acquisition with Adaptive Antenna Array and Adaptive Threshold for DS-CDMA Wireless Communications”, IEEE GLOBECOM, San Francisco, Calif., Nov. 27-Dec. 1, 2000. However, to the inventors' knowledge, practical PNCTLs using all smart antenna elements are not available yet. The invention employs all smart antenna elements for PN code tracking. Verification that tracking performance can be improved significantly when all smart antenna element outputs are added and fed into a typical noncoherent digital PNCTL is set forth herein.

Root-mean-square-chip-time-jitter (RMS) and mean-time-to-lose-lock (MTLL) are typical performance criteria for a PNCTL. Once the PN code is locked, longer MTLL is better for mobile users. Longer MTLL can be achieved as RMS gets smaller. Much analysis and design have appeared for analog and digital phase lock loop (DPLL) and PNCTLs, (e.g., Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995; H. Meyr and G. Ascheid, “Synch. in Digital Communications”, John Wiley and Sons, 1990; R. De Gaudenzi, M. Luise, and R. Viola, “A Digital Chip Timing Recovery Loop for Band Limited Direct-Sequence Spread-Spectrum Signals”, IEEE Trans. on Comm., Vol. 41, no. 11, pp. 1760-1769. November 1993; R. De Gaudenzi and M. Luise, “Decision Directed Coherent Delay Lock Tracking Loop for DS-Spread Spectrum Signals”, IEEE Trans. on Comm., Vol. 39, No. 5, pp. 758-765, May 1991; R. De Gaudenzi, “Direct Sequence Spread Spectrum Chip Tracking in the Presence of Unresolvable Multipath Components”, IEEE Trans. on Vehicular Technology, Vol. 48, No. 5, pp. 1573-1583, September 1999; and A. Weinberg and B. Liu, “Discrete Time Analysis of Nonuniform Sampling First and Second Order Digital Phase Lock Loops”, IEEE Trans. on Comm., Vol. Com-22, No. 2, pp. 123-137, February 1974), under additive white Gaussian noise (AWGN) and fading environments. The analyses in these studies suggest that decreasing loop bandwidth or increasing desired signal transmitting power is the way to improve MTLL. However, the present invention, without changing loop bandwidth and desired signal power, improves MTLL by using all smart antenna elements since a smart antenna can improve the SINR output.

A typical noncoherent early-late gate PNCTL of the first and second order are considered herein (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995.) Then, the analysis in R. De Gaudenzi, M. Luise, and R. Viola, “A Digital Chip Timing Recovery Loop for Band Limited Direct-Sequence Spread-Spectrum Signals”, IEEE Trans. on Comm., Vol. 41, No. 11, pp. 1760-1769, November 1993 and R. De Gaudenzi and M. Luise, “Decision Directed Coherent Delay Lock Tracking Loop for DS-Spread Spectrum Signals”, IEEE Trans. on Comm, Vol. 39, No. 5, pp. 758-765, May 1991, to the second order early-late gate PNCTL, while only the first order was studied in the first of these 2 references. Also, the analysis is extended (A. Weinberg and B. Liu, “Discrete Time Analysis of Nonuniform Sampling First and Second Order Digital Phase Lock Loops”, IEEE Trans. on Comm., Vol. Com-22, No. 2, pp. 123-137, February 1974) to a fading environment and PNCTL because a typical mobile channel is under fading while Weinberg and Liu (in A. Weinberg and B. Liu, “Discrete Time Analysis of Nonuniform Sampling First and Second Order Digital Phase Lock Loops”, IEEE Trans. on Comm., Vol. Com-22, No. 2, pp. 123-137, February 1974) presented the steady state phase error probability density function and variance of the first and second order DPLL only under an AWGN environment.

The proposed PNCTL of smart antennas in a 3G W-CDMA uplink dedicated physical control channel (DPCCH) has been simulated. To do that, the on-time array output was used as input to NLMS adaptive algorithm whose reference value is set to the number of elements, M. The NLMS algorithm requires only (3M+1) complex multiplication and 2M additions and does not require any phase calibration. Array outputs with early and late PN codes are used as inputs to the noncoherent PNCTL. In addition, loop bandwidth and RMS are expressed in terms of poles of the overall PNCTL transfer function for systematic design instead of using a trial and error method. Optimum pole locations to minimize RMS for the first and second order PNCTLs are found. Finally, through both analysis and simulation, significant improvements in RMS and MTLL are demonstrated.

BRIEF DESCRIPTION OF THE DRAWINGS

A more complete appreciation of the invention and many of the attendant advantages and features thereof will be readily understood by reference to the following detailed description when considered in connection with the accompanying drawings wherein:

FIG. 1 is a block diagram of a noncoherent PN code chip time tracking loop (PNCTL) with a smart antenna according to an embodiment of the present invention.

FIG. 2 is an equivalent PN code chip time tracking loop (PNCTL) using a linear model.

FIG. 3 is a plot of the magnitude |{tilde over (f)}(p₁,p₂)| in e.g. (17) versus pole locations.

FIG. 4 is a plot of a normalized loop bandwidth B_(L)T_(c) for the first and second order PN code tracking loops (PNCTLs) using equations (25) and (15), when chip rate is 3.84 Mcps and window size N_(t) is 256 chips for energy accumulation, as a function of pole location.

FIG. 5 is a plot of a root mean square (RMS) chip time jitter of the first order PN code tracking loop under fading and AWGN environments with the number of antenna elements M as a parameter.

FIG. 6 is a plot a root mean square (RMS) chip time jitter of the second order PN code tracking loop under fading and AWGN environments with the number of antenna elements M as a parameter.

FIG. 7 is a plot of the mean time to lose lock (MTLL) of a first order PN code tracking loop under fading and AWGN environments with the number of antenna elements M as a parameter.

FIG. 8 is a plot of the mean time to lose lock (MTLL) of the second order PN code tracking loop under fading and AWGN environments with the number of antenna elements M as a parameter.

DETAILED DESCRIPTION 1. Introduction

The present invention will be discussed with reference to preferred embodiments of methods and devices for psuedonoise code tracking loops. Specific details, such as numbers of antenna elements and the order of loop filters, are set forth in order to provide a through understanding of the present invention. The preferred embodiments discussed herein should not be understood to limit the invention. Furthermore, for ease of understanding, certain method steps are delineated as separate steps; however, these steps should not be construed as necessarily distinct nor order dependent in their performance.

Section 2 describes a noncoherent PNCTL. Section 3 presents analysis: (1) the overall transfer function of PNCTL, and (2) RMS and MTLL of the first and second orders under AWGN and fading environments. Section 4 shows both analysis and simulation results for RMS and MTLL. Section 5 makes conclusions.

2. System Modeling

FIG. 1 shows a block diagram of the proposed noncoherent PNCTL with smart antennas. A single dwell serial PN code search is considered for demonstration. Other search can be employed. When the received and local PN codes are synchronized within a fraction of a chip time, i.e., when PN code acquisition is achieved, the PNCTL in FIG. 1 is triggered. The DPCCH in the W-CDMA, or called 3GPP system (Third Generation Partnership Project, “Spreading and Modulation (FDD)”, 3GPP Technical Specification, TS25.213, v3.2.0, March 2000), is used for the PN code chip time tracking demonstration. The received signal is down-converted and passed through matched filters. In the 3G W-CDMA system, a square root raised cosine (RRC) filter with frequency response H(f) is employed at both the transmitter and receiver with roll-off factor β=0.22. The impulse time response of the filter response H(f) can be written as $\begin{matrix} {{{RRC}(t)} = \frac{{\sin\left( {\pi\frac{t}{T_{c}}\left( {1 - \beta} \right)} \right)} + {4\beta\frac{t}{T_{c}}{\cos\left( {\pi\frac{t}{T_{c}}\left( {1 - \beta} \right)} \right)}}}{\pi\frac{t}{T_{c}}\left( {1 - \left( {4\beta\frac{t}{T_{c}}} \right)^{2}} \right)}} & (1) \end{matrix}$ where T_(c) is a chip time interval. The received signal is multiplied by three local PN codes. One code is advanced by Δ=T_(c)/2 seconds for the early gate, another is delayed by Δ=T_(c)/2 seconds for the late gate, and the third is on-time and used for signal demodulation. The outputs multiplied by the local PN codes are accumulated over N_(t) chips to achieve processing gain, and denoted by an M×1 vector y(k) for M element antenna where k denotes the update index, i.e., kN_(t)T_(c). An inner product y^(H)(k)w(k) is taken at every k to obtain an array output where superscripts H and T in the paper denote the Hermitian and transpose, respectively. The same weight vector w(k) is multiplied to the three de-spread outputs to get array outputs. The on-time array output is fed into an adaptive algorithm to update the weight vector. The NLMS adaptive algorithm is used for weight vector updating (Simon Haykin, “Adaptive Filter Theory”, pp. 437, Prentice Hall 1996). The reference value is set to the number of antenna elements M because ideal weight vector w(k) should be matched with array response vector a(θ(k))=(1, e^(−jπsin(θ(k))), . . . , e^(−j(M−1)π sin(θ(k))))^(T) where θ(k) is the DOA from the desired signal, and ideal array output w^(H)(k)a(θ(k)) would be equal to M. Array outputs for early and late gates are fed into the square devices in the PNCTL, FIG. 1, to measure energy. Measured energy is denoted by Z⁻ and Z⁺ for the advanced and delayed branches, respectively. Energy difference Z_(Δ) is passed through loop filter F(z) and multiplied with loop filter gain γ. The voltage controlled oscillator (VCO) takes the loop filter output and controls the chip time shifting amount for the local PN code. The transfer function of VCO is modeled by z⁻¹/(1−z⁻¹) (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995), where z denotes a sample delay. Four samples are taken per chip. The transfer function of loop filter is denoted by F(z) and ${F(z)} = {1 + \frac{{az}^{- 1}}{1 - z^{- 1}}}$ for a second order PNCTL (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995, pp. 66). The first order PNCTL is a special case of the second order with a=0, where a is for the Doppler rate compensation. The loop parameters a and γ are determined for given loop bandwidths or pole locations. The Jakes Rayleigh fading model is employed for multipath simulation.

3. Analysis

Inverse Fourier transform of |H(f)|² can be written as $\begin{matrix} {{g(\tau)} = {\frac{\cos\left( {{\pi\beta\tau}/T_{c}} \right)}{1 - \left( {2{{\beta\tau}/T_{c}}} \right)^{2}}\frac{\sin\left( {{\pi\tau}/T_{c}} \right)}{{\pi\tau}/T_{c}}}} & (2) \end{matrix}$ where H(f) is the RRC frequency response of time response in (1). The output of the receiver matched filter in FIG. 1 is proportional to g(τ) because the transmitter also employs a square root raised cosine filter H(f). If samples of g(τ) are taken at integer times T_(c), then there would be no degradation in desired signal power when the PN code is synchronized. Let ε≡τ/T_(c) denote the normalized chip timing error. Then, g(τ) at τ=εT_(c) is $\begin{matrix} {{g\left( {ɛ\quad T_{c}} \right)} = {\frac{\cos({\pi\beta ɛ})}{1 - \left( {2{\beta ɛ}} \right)^{2}}{\frac{\sin({\pi ɛ})}{\pi ɛ}.}}} & (3) \end{matrix}$ The loop error characteristic of early and late gate PNCTL (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995) can be written as $\begin{matrix} {{\eta(ɛ)} = {{g^{2}\left\lbrack {\left( {ɛ - \frac{1}{2}} \right)T_{c}} \right\rbrack} - {g^{2}\left\lbrack {\left( {ɛ + \frac{1}{2}} \right)T_{c}} \right\rbrack}}} & (4) \end{matrix}$ by using (3) where early and late gates employ a half chip delayed and advanced version of received samples. The PNCTL characteristic can be assumed to be linear when chip time error ε is small. The slope at ε=0 is obtained by taking a derivative of (4) as $\begin{matrix} \begin{matrix} {{\kappa \equiv \frac{\mathbb{d}\quad{\eta(ɛ)}}{\mathbb{d}\quad ɛ}}❘_{ɛ = 0}} \\ {= {16\quad{\cos\left( \frac{\beta\pi}{2} \right)}\frac{{{{\beta\pi}\left( {1 - \beta^{2}} \right)}{\sin\left( {{\pi\beta}/2} \right)}} + {2\left( {1 - {3\beta^{2}}} \right){\cos\left( {{\pi\beta}/2} \right)}}}{{\pi^{2}\left( {1 - \beta^{2}} \right)}^{3}}}} \\ {\cong 3.2416} \end{matrix} & (5) \end{matrix}$ where β=0.22 is used. Let Z_(Δ) denote the energy difference between the early and late branches in FIG. 1. Then, average energy difference can be written as E[Z _(Δ) ]=N _(t) ² E _(c)η(ε)≅N _(t) ² E _(c)κ  (6) by using the loop error characteristic in (4) (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995).

FIG. 2 shows an equivalent PNCTL model using (6) where V₀ is an upper bound of total interference variance in the tracking loop when a smart antenna is employed and can be written as V ₀=2N _(t) ²(I _(0,in) /M)²(1+2N_(t) E _(c)/(I _(0,in) /M))  (7) where I_(0,in)/2 is variance of input interference fed into each antenna element (Y. S. Song, H. M. Kwon, and B. J. Min, “Computationally Efficient Smart Antennas for CDMA Wireless Comm.”, IEEE Trans. On Vehicular Technology, Vol. 50, No. 6, November 2001; Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995; Yoo S. Song, Hyuck M. Kwon, and Byung J. Min, “Smart Antennas for 3G and Future Generation CDMA”, IEEE International Conference on Phased Array Systems and Technology, Dana Point, Calif., pp. 69-72, May 21-25, 2000). The overall loop transfer function can be found from the equivalent PNCTL as $\begin{matrix} {{H(z)} = {\frac{N_{t}^{2}E_{c}{\kappa\gamma}\quad{F(z)}\frac{z^{- 1}}{1 - z^{- 1}}}{1 + {N_{t}^{2}E_{c}{\kappa\gamma}\quad{F(z)}\frac{z^{- 1}}{1 - z^{- 1}}}} = \frac{N_{t}^{2}E_{c}{\kappa\gamma}\quad{F(z)}z^{- 1}}{1 - z^{- 1} + {N_{t}^{2}E_{c}{\kappa\gamma}\quad{F(z)}z^{- 1}}}}} & (8) \end{matrix}$ which is equal to that of the single antenna element case (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, Addison Wesley, N.Y., 1995). The major difference between a smart antenna and a single antenna element is interference suppressed by a factor M as (7). 3.1. Root Mean Square Chip Time Jitter (RMS) Analysis

By using ${F(z)} = {1 + \frac{{az}^{- 1}}{1 - z^{- 1}}}$ for the second order PNCTL, (8) can be rewritten as $\begin{matrix} {{H(z)} = {\frac{N_{t}^{2}E_{c}{\kappa\gamma}\quad\left( {z^{- 1} - {\left( {1 - a} \right)z^{- 2}}} \right)}{1 + {\left( {{- 2} + {N_{t}^{2}E_{c}{\kappa\gamma}}}\quad \right)z^{- 1}} + {\left( {1 - {N_{t}^{2}E_{c}{\kappa\gamma}\quad\left( {1 - a} \right)}} \right)z^{- 2}}}.}} & (9) \end{matrix}$ Let p₁ and p₂ denote two poles of H(z) in (9). Then, the two poles satisfy the following: p ₁ +p ₂=2−N _(t) ² E _(c)κγ  (10-a) and p ₁ p ₂=1−N _(t) ² E _(c)κγ(1−a).  (10-b) Therefore, loop filter gain γ and Doppler rate compensation coefficient a can be expressed in terms of poles as $\begin{matrix} {{\gamma = \frac{2 - p_{1} - p_{2}}{N_{t}^{2}E_{c}\kappa}}{and}} & \left( {11\text{-}a} \right) \\ {a = {\frac{1 - p_{1} - p_{1} + {p_{1}p_{2}}}{2 - p_{1} - p_{2}}.}} & \left( {11\text{-}b} \right) \end{matrix}$ Variance of chip time error for the second order PNCTL is derived in Appendix in detail and can be written as $\begin{matrix} \begin{matrix} {{{Var}\left( \frac{\tau}{T_{c}} \right)} = {\frac{V_{0}}{\left( {N_{t}^{2}E_{c}\kappa} \right)^{2}}{\oint{{H(z)}{H\left( z^{- 1} \right)}\frac{\mathbb{d}z}{2{\pi\mathbb{i}z}}}}}} \\ {= {\frac{V_{0}\gamma^{2}}{\left( {N_{t}^{2}E_{c}{\kappa\gamma}} \right)^{2}\left( {a - 1} \right)}\frac{{2a} + {2N_{t}^{2}E_{c}{\kappa\gamma}} - {3{aN}_{t}^{2}E_{c}{\kappa\gamma}} + {a^{2}N_{t}^{2}E_{c}{\kappa\gamma}}}{\left( {{2N_{t}^{2}E_{c}{\kappa\gamma}} - {{aN}_{t}^{2}E_{c}{\kappa\gamma}} - 4} \right)}}} \end{matrix} & (12) \end{matrix}$ while variance of the first order PNCTL is shown in (Andrew J. Viterbi, “CDMA Principle of Spread Spectrum Communication”, pp. 66, Addison Wesley, N.Y., 1995) which is a special case of (12) with a=0. One-side loop bandwidth in Hertz is given as $\begin{matrix} {B_{L} = {\int_{0}^{\pi/T_{c}}{\left\lbrack {{H\left( {\mathbb{e}}^{{j\omega}\quad T_{c}} \right)}{H\left( {\mathbb{e}}^{{- {j\omega}}\quad T_{c}} \right)}} \right\rbrack{\frac{\quad{\mathbb{d}\omega}}{2\pi}.}}}} & (13) \end{matrix}$ Thus, (12) can be expressed in terms of loop bandwidth by using (13) as $\begin{matrix} {{{Var}\left( \frac{\tau}{T_{c}} \right)} = {\frac{2V_{0}B_{L}T_{c}}{\left( {N_{t}^{2}E_{c}\kappa} \right)^{2}}.}} & (14) \end{matrix}$ Normalized loop bandwidth of the second order PNCTL can be obtained from (12) and (14) and expressed in terms of poles with (10) as $\begin{matrix} \begin{matrix} {{B_{L}T_{c}} = \frac{{2a} + {2N_{t}^{2}E_{c}{\kappa\gamma}} - {3{aN}_{t}^{2}E_{c}{\kappa\gamma}} + {a^{2}N_{t}^{2}E_{c}{\kappa\gamma}}}{2\left( {a - 1} \right)\left( {{2N_{t}^{2}E_{c}{\kappa\gamma}} - {{aN}_{t}^{2}E_{c}{\kappa\gamma}} - 4} \right)}} \\ {= {\frac{{2\left( {2 - p_{1} - p_{2}} \right)} + {a\left( {{p_{1}p_{2}} + {2p_{1}} + {2p_{2}} - 3} \right)}}{{- 2}\left( {a - 1} \right)\left( {{p_{1}p_{2}} + p_{1} + p_{2} + 1} \right)}.}} \end{matrix} & (15) \end{matrix}$ RMS of the second order PNCTL can be obtained by taking the square root of (12) or (14). However, it is more useful to represent RMS in terms of loop bandwidth as $\begin{matrix} {\sigma_{ɛ} = {\frac{\sqrt{2V_{0}B_{L}T_{c}}}{N_{t}^{2}E_{c}\kappa}.}} & (16) \end{matrix}$ Equation (16) shows that RMS is proportional to loop bandwidth B_(L)T_(c).

To find the optimum poles that minimize variance of chip time error, we take a derivative of (12) with respect to p₁ and p₂, after substituting (11) into (12) as $\begin{matrix} {\begin{matrix} {{{\frac{\partial^{2}}{{\partial p_{1}}{\partial p_{2}}}{{Var}\left( \frac{\tau}{T_{c}} \right)}}❘_{{\gamma = \frac{2 - p_{1} - p_{2}}{N_{t}^{2}E_{c}\kappa}},{a = \frac{1 - p_{1} - p_{2} + {p_{1}p_{2}}}{2 - p_{1} - p_{2}}}}} = {\frac{V_{0}}{\left( {N_{t}^{2}E_{c}\kappa} \right)^{2}}{\overset{\sim}{f}\left( {p_{1},p_{2}} \right)}}} \\ {= 0} \end{matrix}{where}} & (17) \\ {{\overset{\sim}{f}\left( {p_{1},p_{2}} \right)} = \frac{{Num}_{\overset{\sim}{f}}}{{Den}_{\overset{\sim}{f}}}} & (18) \end{matrix}$ Num _({tilde over (f)}) =p ₁ ³ p ₂ ³+2p ₁ ² p ₂ ³ +p ₁ p ₂ ³+2p ₁ p ₂ ³+2p ₁ ³ p ₂ ²−11p ₁ ² p ₂ ²−4p ₁ p ₂ ² +p ₂ ² +p ₂ ² +p ₁ ³ p ₂−4p ₁ ² p ₂+13p ₁ p ₂ +p ₁ ²+2p ₁−7  (19) Den _({tilde over (f)})=0.5(1+p ₁)(1+p ₂)(p ₁ p ₂−1)³(1+p ₁ +p ₂ +p ₁ p ₂).  (20) FIG. 3 shows magnitude |{tilde over (f)}(p₁,p₂)| by varying two pole locations, p₁ and p₂. If two poles p₁ and p₂ are located in the right hand side of unit circle z-domain, it is observed that |{tilde over (f)}(p₁,p₂)| approaches zero as poles are close to 1 on the real axis. This implies that optimum pole locations are equal to 1. Therefore, variance of chip time error will be minimum as poles approach 1 on the real axis. However, pole locations are not allowed to exceed, nor to equal, 1 because of tracking loop stability.

Under fading environment, (7) can be modified into conditional variance of interference V_(0|α) for given fading amplitude a as V _(0|α)2N _(t) ²(I _(0,in) /M)²+4α² N _(t) ³(I _(0,in) /M)E _(c).  (21) Then, the average variance of chip timing error under fading environment can be obtained with a probability density function (pdf) of fading amplitude, p_(R)(α), as $\begin{matrix} {{{{Var}\left( \frac{\tau}{T_{c}} \right)} = {\int_{0}^{\infty}{{{Var}\left( {\frac{\tau}{T_{c}}❘\alpha} \right)}{p_{R}(\alpha)}\quad{\mathbb{d}\alpha}}}}{where}} & (22) \\ \begin{matrix} {{{Var}\left( {\frac{\tau}{T_{c}}❘\alpha} \right)} = {\frac{V_{0|\alpha}}{\left( {N_{t}^{2}\alpha^{2}E_{c}} \right)^{2}\kappa^{2}}{\oint{{H\left( z \middle| \alpha \right)}{H\left( z^{- 1} \middle| \alpha \right)}\frac{\mathbb{d}z}{2\pi\quad{\mathbb{i}z}}}}}} \\ {= \frac{V_{0|\alpha}}{\left( {N_{t}^{2}\alpha^{2}E_{c}\kappa} \right)^{2}\left( {a - 1} \right)}} \\ {\frac{{2a} + {2N_{t}^{2}\alpha^{2}E_{c}{\kappa\gamma}} - {3{aN}_{t}^{2}\alpha^{2}E_{c}{\kappa\gamma}} + {a^{2}N_{t}^{2}\alpha^{2}E_{c}{\kappa\gamma}}}{{2N_{t}^{2}\alpha^{2}E_{c}{\kappa\gamma}} - {{aN}_{t}^{2}\alpha^{2}E_{c}{\kappa\gamma}} - 4}} \end{matrix} & (23) \end{matrix}$ For a Rayleigh fading model, $\begin{matrix} {{p_{R}(\alpha)} = {\frac{\alpha}{\sigma_{\alpha}^{2}}{\mathbb{e}}^{\frac{- \alpha^{2}}{2\sigma_{\alpha}^{2}}}}} & (24) \end{matrix}$ where σ_(α) ²=1. RMS under fading channel can be obtained by taking a square root of (22).

Similarly, analysis for the first order PNCTL can be done as a special case of the second order with a=0 in (9)-(23). The first order PNCTL transfer function has a single pole p at 1−N_(t) ²E_(c)κγ, and loop filter gain γ can be expressed in terms of the pole as γ=(1−p)/N_(t) ²E_(c)κ. And normalized loop bandwidth of the first order PNCTL can be written as $\begin{matrix} {{B_{L}T_{c}} = {\frac{N_{t}^{2}E_{c}{\kappa\gamma}}{2\left( {2 - {N_{t}^{2}E_{c}{\kappa\gamma}}} \right)} = {\frac{1 - p}{2\left( {1 + p} \right)}.}}} & (25) \end{matrix}$ RMS of the first order PNCTL can be obtained by replacing B_(L)T_(c) in (14) with (25) and taking a square root of (14). The optimum pole location that minimizes variance (or RMS) can be found by replacing γ with (1−p)/N_(t) ²E_(c)κ and taking a derivative of variance with respect to p as $\begin{matrix} {{\frac{\mathbb{d}}{\mathbb{d}p}\left( \left. \frac{V_{0}\gamma}{N_{t}^{2}E_{c}{\kappa\left( {2 - {N_{t}^{2}E_{c}{\kappa\gamma}}} \right)}} \right|_{\gamma = \frac{1 - p}{N_{t}^{2}E_{c}\kappa}} \right)} = {\frac{V_{0}\left( {- 2} \right)}{\left( {N_{t}^{2}E_{c}\kappa} \right)^{2}\left( {1 + p} \right)^{2}} < 0.}} & (26) \end{matrix}$ Note from (26) that the derivative is negative, which implies that variance decreases asp increases. Therefore, RMS will be minimum when p is 1. RMS of the first order PNCTL under fading environment can be obtained similar to (22). 3.2. Mean Time to Lose Lock (MTLL) Analysis

Loop noise Ar(k) of variance V₀ in FIG. 2 is the difference between Z_(Δ)(k) and its average E[Z_(Δ)(k)] and can be written as N(k)=Z _(Δ)(k)−E[Z _(Δ)(k)]=Z _(Δ)(k)−N _(t) ² E _(c)η{ε(k)}  (27) from (6) where k denotes the PNCTL update index. Since loop filter F(z) can be represented as 1+az⁻¹/(1-z⁻¹) for the second order PNCTL and ε(k) is the VCO output in FIG. 1, the z-transform of ε(k) can be written as $\begin{matrix} {{ɛ(z)} = {{\gamma\left( {- \frac{z^{- 1}}{1 - z^{- 1}}} \right)}\left( {1 + \frac{a\quad z^{- 1}}{1 - z^{- 1}}} \right){Z_{\Delta}(z)}}} & (28) \end{matrix}$ where the negative sign is used in (28) because the VCO output is subtracted from the input timing in FIG. 2. Z_(Δ)(z) in (28) denotes the forward z-transform for Z_(Δ)(k). By taking the inverse z-transform, (28) can be rewritten as $\begin{matrix} \begin{matrix} {{{ɛ(k)} - {2{ɛ\left( {k - 1} \right)}} + {ɛ\left( {k - 2} \right)}} = {- {\gamma\left\lbrack {{Z_{\Delta}\left( {k - 1} \right)} + {\left( {a - 1} \right){Z_{\Delta}\left( {k - 2} \right)}}} \right\rbrack}}} \\ {= {{{- \gamma}\quad N_{\iota}^{2}E_{c}\eta\left\{ {ɛ\left( {k - 1} \right)} \right\}} -}} \\ {{\gamma\quad N\left( {k - 1} \right)} - \left( {a - 1} \right)} \\ {\left\lbrack {{\gamma\quad N_{\iota}^{2}E_{c}\eta\left\{ {ɛ\left( {k - 2} \right)} \right\}} + {\gamma\quad N\left( {k - 2} \right)}} \right\rbrack.} \end{matrix} & (29) \end{matrix}$ Let ε(k)≡s(k+1)+(a−1)s(k). Then (29) becomes $\begin{matrix} {{\left\lbrack {{s\left( {k + 1} \right)} - {2{s(k)}} + {s\left( {k - 1} \right)}} \right\rbrack + {\left( {a - 1} \right)\left\lbrack {{s(k)} - {2{s\left( {k - 1} \right)}} + {s\left( {k - 2} \right)}} \right\rbrack}} = \quad{{{- \gamma}\quad N_{\iota}^{2}E_{c}\eta\left\{ {{s(k)} + {\left( {a - 1} \right){s\left( {k - 1} \right)}}} \right\}} - {\gamma\quad{N\left( {k - 1} \right)}} + {{\left( {a - 1} \right)\left\lbrack {{{- \gamma}\quad N_{\iota}^{2}E_{c}\eta\left\{ {{s\left( {k - 1} \right)} + {\left( {a - 1} \right){s\left( {k - 2} \right)}}} \right\}} - {\gamma\quad{N\left( {k - 2} \right)}}} \right\rbrack}.}}} & (30) \end{matrix}$ Equation (30j can be broken into two equivalent equations. Let y₁(k)=s(k−2) and y₂(k)=s(k−1). Then, (30) can be rewritten as the following state equations: y ₁(k+1)=y ₂(k) y ₂(k+1)=2y ₂(k)−y ₁(k)−γN_(t) ² E _(c)η(y ₂(k)+(a−1)y ₁(k))−γN(k−2)  (31) where ε(k)=y ₂(k+2)+(a−1)y ₁(k+2).  (32) The y₁(k) and y₂(k) are discrete time continuous variable Markov processes. The joint pdf of y₁(k) and y₂(k) satisfies the Chapman-Kolmogorov equation (A. Weinberg and B. Liu, “Discrete Time Analysis of Nonuniform Sampling First and Second Order Digital Phase Lock Loops,” IEEE Trans. on Comm., Vol. Com-22, No. 2, pp. 123-137, February 1974) as $\begin{matrix} {{p_{k + 1}\left( {y_{1},\left. y_{2} \middle| y_{1,0} \right.,y_{2,0}} \right)} = {\int_{\quad^{- \infty}}^{\infty}{\int_{- \infty}^{\infty}{{q_{k}\left( {y_{1},\left. y_{2} \middle| x_{1} \right.,x_{2}} \right)}{p_{k}\left( {x_{1},\left. x_{2} \middle| y_{1,0} \right.,y_{2,0}} \right)}{\mathbb{d}x_{1}}{\mathbb{d}x_{2}}}}}} & (33) \end{matrix}$ where y_(1,0)=y₁(0) and y_(2,0)=y₂(0) are the initial chip timing errors, p_(k)(·,·|y_(1,0), y_(2,0)) is the joint pdf of y₁(k) and y₂(k) for given y_(1,0) and y_(2,0), and q_(k)(y₁, y₂|x₁, x₂) is the joint transition pdf of y₁(k+1) and y₂(k+1) for given y₁(k)=x₁ and y₂(k)=x₂. Noting that noise N(k) is independent of y₁(k) and y₂(k), we can observe from (31) that y₂(k+1) for given (y₁(k)=x₁, y₂(k)=x₂) is Gaussian with mean and variance equal to E_(k)(y₁, y₂|x₁, x₂)=2x₂−x₁−γN_(t) ²E_(c)η(x₂+(a−1)x₁) and γ²V₀, respectively. Thus, the conditional density of y₂(k+1), given (y₁(k)=x₁, y₂(k)=x₂), is independent of k since the channel is assumed to be stationary. Therefore, from (31), one obtains $\begin{matrix} {{q_{k}\left( {y_{1},\left. y_{2} \middle| x_{1} \right.,x_{2}} \right)} = {{\delta\left( {y_{1} - x_{2}} \right)}\frac{1}{\sqrt{2{\pi\gamma}^{2}V_{0}}}{{\exp\left\lbrack {- \frac{\left( {y_{2} - {2x_{2}} + x_{1} + {\gamma\quad N_{\iota}^{2}E_{c}{\eta\left( {x_{2} + {\left( {a - 1} \right)x_{1}}} \right)}}} \right)^{2}}{2\gamma^{2}V_{0}}} \right\rbrack}.}}} & (34) \end{matrix}$ Then (33) becomes $\begin{matrix} {\begin{matrix} {{p_{k + 1}\left( {y_{1},y_{2}} \right)} = {\int_{- \infty}^{\infty}{\frac{1}{\sqrt{2{\pi\gamma}^{2}V_{0}}}\exp}}} \\ {\left\lbrack {- \frac{\left( {y_{2} - {2y_{1}} + x_{1} + {\gamma\quad N_{\iota}^{2}E_{c}{\eta\left( {y_{1} + {\left( {a - 1} \right)x_{1}}} \right)}}} \right)^{2}}{2\gamma^{2}V_{0}}} \right\rbrack} \\ {{p_{k}\left( {x_{1},y_{1}} \right)}{\mathbb{d}x_{1}}} \\ {= {\int_{- \infty}^{\infty}{{K\left( {y_{1},y_{2},x_{1}} \right)}{p_{k}\left( {x_{1},y_{1}} \right)}{\mathbb{d}x_{1}}}}} \end{matrix}{where}} & (35) \\ {{K\quad\left( {y_{\quad 1},y_{\quad 2},x_{\quad 1}} \right)} = {\frac{1}{\sqrt{2{\pi\gamma}^{2}V_{0}}}{{\exp\left\lbrack {- \frac{\left( {y_{2} - {2y_{1}} + x_{1} + {\gamma\quad N_{\iota}^{2}E_{c}{\eta\left( {y_{1} + {\left( {a - 1} \right)x_{1}}} \right)}}} \right)^{2}}{2\gamma^{2}V_{0}}} \right\rbrack}.}}} & (36) \end{matrix}$ From (31), and since N(k−1) is Gaussian with zero mean and variance V₀, y₁(k+1) and y₂(k+1) will be also jointly Gaussian for all k with zero means. The jointly Gaussian density function is written as $\begin{matrix} {{p_{k + 1}\left( {y_{1},y_{2}} \right)} = {\frac{1}{2{{\pi\sigma}^{2}\left( {k + 1} \right)}\sqrt{1 - {\rho^{2}\left( {k + 1} \right)}}}\exp\left\{ {- \frac{y_{1}^{2} - {2{\rho\left( {k + 1} \right)}y_{1}y_{2}} + y_{2}^{2}}{{2\left\lbrack {1 - {\rho^{2}\left( {k + 1} \right)}} \right\rbrack}{\sigma^{2}\left( {k + 1} \right)}}} \right\}}} & (37) \end{matrix}$ where σ²(k+1) and ρ(k+1) denote variances of y₁(k+1) and y₂(k+1) and their correlation coefficients, respectively. The ρ(k+1) and σ(k+1) are two unknown variables and can be obtained from (35) and (37) by computing two pair points (y₁, y₂) recursively for each k and assuming that the digital loop is initially locked, i.e., y_(1,0)=0 and y_(2,0)=0, i.e., p₀(y₁,y₂)=δ(y₁,y₂), where δ(x) is the Dirac delta function. For example, $\begin{matrix} {\begin{matrix} {{p_{k + 1}\left( {0,0} \right)} = \frac{1}{2{{\pi\sigma}^{2}\left( {k + 1} \right)}\sqrt{1 - {\rho^{2}\left( {k + 1} \right)}}}} \\ {= {\frac{1}{2{{\pi\sigma}^{2}(k)}\sqrt{2{\pi\gamma}^{2}V_{0}}\sqrt{1 - {\rho^{2}(k)}}}{\int_{- \infty}^{\infty}\exp}}} \\ {\left\lbrack {{- \frac{\left( {x_{1} + {\gamma\quad N_{\iota}^{2}E_{c}\eta\left\{ {\left( {a - 1} \right)x_{1}} \right\}}} \right)^{2}}{2\gamma^{2}V_{0}}} - \frac{x_{1}^{2}}{{2\left\lbrack {1 - {\rho^{2}(k)}} \right\rbrack}{\sigma^{2}(k)}}} \right\rbrack{\mathbb{d}x_{1}}} \end{matrix}{and}} & (38) \\ \begin{matrix} {{p_{k + 1}\left( {n,0} \right)} = {\frac{1}{2{{\pi\sigma}^{2}\left( {k + 1} \right)}\sqrt{1 - {\rho^{2}\left( {k + 1} \right)}}}\exp}} \\ {\left\{ {- \frac{n^{2}}{{2\left\lbrack {1 - {\rho^{2}\left( {k + 1} \right)}} \right\rbrack}{\sigma^{2}\left( {k + 1} \right)}}} \right\}} \\ {= {\frac{1}{2{{\pi\sigma}^{2}(k)}\sqrt{2{\pi\gamma}^{2}V_{0}}\sqrt{1 - {\rho^{2}(k)}}}{\int_{- \infty}^{\infty}\exp}}} \\ {\begin{bmatrix} {{- \frac{\left( {{{- 2}n} + x_{1} + {\gamma\quad N_{\iota}^{2}E_{c}{\eta\left( {n + {\left( {a - 1} \right)x_{1}}} \right)}}} \right)^{2}}{2\gamma^{2}V_{0}}} -} \\ \frac{x_{1}^{2} - {2{\rho(k)}x_{1}n} + n^{2}}{2\left( {1 - {\rho^{2}(k)}} \right){\sigma^{2}(k)}} \end{bmatrix}{\mathbb{d}x_{1}}} \end{matrix} & (39) \end{matrix}$ where n is any nonzero real number and set to 0.01. We compute the rightmost terms in (38) and (39) numerically and recursively for each k to obtain p_(k+1)(0,0) and p_(k+1)(n,0). Thus, we have $\begin{matrix} {{\sigma^{2}\left( {k + 1} \right)} = \frac{1}{2\pi\sqrt{1 - {\rho^{2}\left( {k + 1} \right)}}{p_{\quad{k\quad + \quad 1}}\left( {0,0} \right)}}} & (40) \\ {{p_{k + 1}\left( {n,0} \right)} = {{p_{k + 1}\left( {0,0} \right)}\exp\left\{ {- \frac{n^{2}\pi\quad{p_{k + 1}\left( {0,0} \right)}}{\sqrt{1 - {\rho^{2}\left( {k + 1} \right)}}}} \right\}}} & (41) \end{matrix}$ Using (40) and (41), we can solve ρ(k+1) as $\begin{matrix} {{{\rho\left( {k + 1} \right)} = \frac{\sqrt{{- 4}{AC}}}{2A}},{{{\rho\left( {k + 1} \right)}} < 1.}} & (42) \end{matrix}$ where ${A = {\left\lbrack {\ln\left( \frac{p_{k}\left( {n,0} \right)}{p_{k + 1}\left( {0,0} \right)} \right)} \right\rbrack^{2}\quad{and}}}\quad$ $C = {{\left( n^{2} \right)^{2}\pi^{2}{p_{k + 1}^{2}\left( {0,0} \right)}} - {\left\lbrack {\ln\left( \frac{p_{k}\left( {n,0} \right)}{p_{k + 1}\left( {0,0} \right)} \right)} \right\rbrack^{2}.}}$ Also, we can have σ²(k+1) by substituting ρ(k+1)in (42) into (40).

The mean and variance of timing error ε(k) can be obtained from (32) as $\begin{matrix} {{\mu_{ɛ,k} = {{\int_{- \infty}^{\infty}{\int_{- \infty}^{\infty}{\left\lbrack {y_{2} + {\left( {a - 1} \right)y_{1}}} \right\rbrack{p_{k + 2}\left( {y_{1},y_{2}} \right)}\quad{\mathbb{d}y_{1}}\quad{\mathbb{d}y_{2}}}}} = 0}}{and}} & (43) \\ \begin{matrix} {\sigma_{ɛ,k}^{2} = {\int_{- \infty}^{\infty}{\int_{- \infty}^{\infty}{\left\lbrack {y_{2} + {\left( {a - 1} \right)y_{1}} - \mu_{ɛ,k}} \right\rbrack^{2}{p_{k + 2}\left( {y_{1},y_{2}} \right)}\quad{\mathbb{d}y_{1}}\quad{\mathbb{d}y_{2}}}}}} \\ {= {{\sigma^{2}\left( {k + 2} \right)} + {2\left( {a - 1} \right){\sigma^{2}\left( {k + 2} \right)}{\rho\left( {k + 2} \right)}} +}} \\ {\left( {a\quad - \quad 1} \right)^{2}{{\sigma^{\quad 2}\left( {k + 2} \right)}.}} \end{matrix} & (44) \end{matrix}$ Therefore, the pdf of chip timing error ε(k) can be written as $\begin{matrix} {{p_{k}(ɛ)} = {\frac{1}{\sqrt{2\quad\pi\quad\sigma_{ɛ,k}^{2}}}{{\exp\left\lbrack {- \frac{\left( {ɛ - \mu_{ɛ,k}} \right)^{2}}{2\quad\sigma_{ɛ,k}^{2}}} \right\rbrack}.}}} & (45) \end{matrix}$ The probability that the timing error is out of the limit±ε_(L) for the first time at the k-th iteration is denoted by P_(k) and can be written as P _(k) =Pr{# of iterations=k for which |ε_(k)|≧ε_(L) }=Pr{# of iterations≦k for which|ε_(k)|≧ε_(L) }−Pr{of iterations<k for which |ε_(k)|≧ε_(L)} $\begin{matrix} \begin{matrix} {P_{k} = {\Pr\left\{ {{\#\quad{of}\quad{iterations}} = {{k\quad{for}\quad{which}\quad{ɛ_{k}}} \geq ɛ_{L}}} \right\}}} \\ {= {{\Pr\left\{ {{\#\quad{of}\quad{iterations}} \leq {k\quad{for}\quad{which}\quad{ɛ_{k}}} \geq ɛ_{L}} \right\}} -}} \\ {\Pr\left\{ {{\#\quad{of}\quad{iterations}} < {k\quad{for}\quad{which}\quad{ɛ_{k}}} \geq ɛ_{L}} \right\}} \\ {= {\left( {1 - {\int_{- ɛ_{L}}^{ɛ_{L}}{{p_{k}(ɛ)}\quad{\mathbb{d}ɛ}}}} \right) - \left( {1 - {\int_{- ɛ_{L}}^{ɛ_{L}}{{p_{k - 1}(ɛ)}\quad{\mathbb{d}ɛ}}}} \right)}} \\ {= {{\int_{- ɛ_{L}}^{ɛ_{L}}{{p_{k - 1}(ɛ)}\quad{\mathbb{d}ɛ}}} - {\int_{- ɛ_{L}}^{ɛ_{L}}{{p_{k}(ɛ)}\quad{{\mathbb{d}ɛ}.}}}}} \end{matrix} & (46) \end{matrix}$ Therefore, MTLL will be $\begin{matrix} \begin{matrix} {{MTLL} = {\sum\limits_{k = 1}^{\infty}{kP}_{k}}} \\ {= {\sum\limits_{k = 1}^{\infty}{k\left( {\int_{- ɛ_{L}}^{ɛ_{L}}{\left\lbrack {{p_{k - 1}(ɛ)} - {p_{k}(ɛ)}}\quad \right\rbrack{\mathbb{d}ɛ}}} \right)}}} \\ {= {1 + {\sum\limits_{k = 1}^{\infty}{\int_{- ɛ_{L}}^{ɛ_{L}}{{p_{k}(ɛ)}\quad{{\mathbb{d}ɛ}.}}}}}} \end{matrix} & (47) \end{matrix}$

Under fading environment, first we obtain the conditional probability density of the timing error ε(k) for a given fading amplitude α similarly to (45), and then take an average of the conditional pdf. Thus, the MTLL can be obtained under fading environment. For example, under fading environment, the steady state pdf of ε can be written as $\begin{matrix} {{{p(ɛ)} = {\int_{0}^{\infty}{{\alpha\left( {\frac{1}{\sqrt{2\quad\pi\quad\sigma_{ɛ{\alpha}}^{2}}}{\exp\left\lbrack {- \frac{ɛ^{2}}{2\quad\sigma_{ɛ{\alpha}}^{2}}} \right\rbrack}} \right)}{\mathbb{e}}^{\frac{- \alpha^{2}}{2}}\quad{\mathbb{d}\alpha}}}}{where}} & (48) \\ {\sigma_{ɛ{\alpha}}^{2} = {\int_{- \infty}^{\infty}{\int_{- \infty}^{\infty}{\left\lbrack {y_{2} + {\left( {a - 1} \right)y_{1}}} \right\rbrack^{2}{p\left( {y_{1},{y_{2}\left. \alpha \right)\quad{\mathbb{d}y_{1}}\quad{\mathbb{d}y_{2}}}} \right.}}}}} & (49) \\ \begin{matrix} {{p\text{(}y_{1}},{{y_{2}\left. \alpha \right)} = {\int_{- \infty}^{\infty}{\frac{1}{\sqrt{2\quad\pi\quad\gamma^{2}V_{0{\alpha}}}}\exp}}}} \\ {\left\lbrack {- \frac{\begin{matrix} \left( {y_{2} - {2y_{1}} + x + {\gamma\quad\alpha^{2}N_{t}^{2}}} \right. \\ \left. {E_{c}{\eta\left( {y_{1} + {\left( {a - 1} \right)x}} \right)}} \right)^{2} \end{matrix}}{2\quad\gamma^{2}V_{0{\alpha}}}} \right\rbrack{p\left( {x,y_{1}} \right)}{\mathbb{d}x}} \end{matrix} & (50) \end{matrix}$ and V_(0|a) is given in (21).

The MTLL for the first order PNCTL in (R. DeGaudenzi, M. Luise, and R. Viola, “A Digital Chip Timing Recovery Loop for Band Limited Direct-Sequence Spread-Spectrum Signals,” IEEE Trans. on Comm., Vol. 41, No. 11, pp. 1760-1769, November 1993) is a special case of the second order PNCTL described herein, with a=0.

4. Numerical and Simulation Results

For simulation, the DOA of the desired signal was chosen arbitrarily between −90° and 90° with a zero angle spread. Interference was added at each antenna element with a mean of zero and a variance of I_(0,in)/2. A linear array of M=4 antenna elements was assumed. The constant parameter and convergence parameter in the NLMS were chosen as 0.1 and 1.5, respectively. The number of chip intervals N_(t) for energy accumulation was set to spreading factor 256 of DPCCH. The antenna weight vector and PNCTL were updated at every control symbol time=0.006 ms. Normalized bandwidth B_(L)T_(c) was chosen as 0.0556. The corresponding pole for the first order PNCTL was 0.8. The corresponding pole locations for the second order PNCTLs were 0.9 and 0.9364. Also the RRC filter specified in W-CDMA was used, and so κ=3.2.

Jakes fading model was used with mobile speeds in the range of 10 km/h and 1000 km/h. Although results are not shown in this paper, it was observed that RMS chip time jitter and MTLL are not sensitive to mobile speeds. This may be due to noncoherent energy detection in the tracking loop.

FIG. 4 shows normalized bandwidth B_(L)T_(c) versus the second pole location¹ with the first pole location as a parameter using (25) and (15) for the second order PNCTL, respectively. As the pole approaches 1, the bandwidth gets smaller, which means RMS decreases according to (16) and MTLL increases as indicated in (R. DeGaudenzi, M. Luise, and R. Viola, “A Digital Chip Timing Recovery Loop for Band Limited Direct-Sequence Spread-Spectrum Signals,” IEEE Trans. on Comm., Vol. 41, No. 11, pp. 1760-1769, November 1993). ¹ For the first order PNCTL, the pole location in the horizontal axis of FIG. 4 represents its single pole location.

FIGS. 5 and 6 show RMS under AWGN and fading environments for the first and second order PNCTL, respectively, with the number of smart antenna elements M as a parameter. The solid and dashed lines represent analysis results under AWGN and fading environments, respectively. The triangles and circles show corresponding simulation results for AWGN and fading environments, respectively. Results show significant improvement by employing smart antennas, e.g., 6 dB improvement in E_(c)/I₀ with M=4 elements to achieve the same RMS, compared to M=1. Performance under fading environment is slightly degraded, compared to AWGN, especially at low E_(c)/I₀. It is also observed that simulation results agree well with analysis.

FIGS. 7 and 8 show MTLL under AWGN and fading environments for the first and second order PNCTLs, respectively, with M as a parameter. Again, performance is significantly improved, e.g., by 6 dB with a smart antenna of M=4 elements. Fading channel degrades performance by 0.5 dB, compared to AWGN. The second order PNCTL is about 1 dB better than the first order for the same bandwidth.

5. Conclusions

The invention employs a smart antenna to improve a noncoherent PN code tracking loop. RMS and MTLL were analyzed for the first and second order PNCTLs under fading and AWGN environments. Simulation results agree well with the theoretical results. The proposed PN tracking scheme was simulated with a 3G W-CDMA DPCCH channel. Chip time tracking performance can be improved significantly, e.g., by 6 dB with four antenna elements, compared to a tracking scheme with a single antenna element. For a given bandwidth, the chip time tracking loop with the second order PNCTL can achieve a better MTLL than that achieves with the first order. RMS and MTLL performance under fading environment are slightly worse than that under AWGN. In addition, relations between poles of the PNCTL transfer function with normalized bandwidth and RMS, allow the use of a systematic PNCTL design instead of the tedious trial and error method. Optimum pole locations are 1 on the real axis to minimize RMS for both the first and second order PNCTLs. In conclusion, a smart antenna can significantly improve chip time tracking performance as well as bit error rate.

Obviously, numerous other modifications and variations of the present invention are possible in light of the above teachings. It is therefore to be understood that within the scope of the appended claims, the invention may be practiced otherwise than as specifically described herein.

APPENDIX

Variance of Chip Timing Error in Second Order PN Code Tracking Loop. $\begin{matrix} {\begin{matrix} {{{Var}\quad\left( \quad\frac{\tau}{\quad T_{\quad c}} \right)} = {\frac{V_{0}}{\left( {N_{t}^{2}E_{c}\kappa} \right)^{2}}{\oint{{H(z)}{H\left( z^{- 1} \right)}\frac{\mathbb{d}z}{2\quad\pi\quad{\mathbb{i}}\quad z}}}}} \\ {= {\frac{V_{0}\gamma^{2}}{2\quad\pi\quad{\mathbb{i}}}{\oint{{q(z)}{\mathbb{d}z}}}}} \end{matrix}{where}} & ({A1}) \\ \begin{matrix} {{q(z)} = \frac{\left\lbrack {z + \left( {{- 1} + a} \right)} \right\rbrack\left\lbrack {1 + {\left( {{- 1} + a} \right)z}} \right\rbrack}{\left\lbrack {z^{2} + {Az} + B} \right\rbrack\left\lbrack {1 + {Az} + {Bz}^{2}} \right\rbrack}} \\ {= \frac{\left\lbrack {z + \left( {{- 1} + a} \right)} \right\rbrack\left\lbrack {1 + {\left( {{- 1} + a} \right)z}} \right\rbrack}{\left( {z - z_{1}} \right)\left( {z - z_{2}} \right)\left( {1 - {z_{1}z}} \right)\left( {1 - {z_{2}z}} \right)}} \end{matrix} & ({A2}) \\ {A = {{- 2} + {N_{t}^{2}E_{c}\kappa\quad\gamma}}} & ({A3}) \\ {B = {1 - {N_{t}^{2}E_{c}\kappa\quad{{\gamma\left( {1 - a} \right)}.}}}} & ({A4}) \end{matrix}$ Let z₁ and z₂ denote roots of z²+Az+B=0 inside the unit circle. Then 1/z₁ and 1/z₂ are roots of 1+Az+Bz²=0, which are outside the unit circle. We know z₁+z₂=−A and z₁z₂=B. Then, $\begin{matrix} \begin{matrix} {{{Var}\quad\left( \frac{\tau}{\quad T_{\quad c}} \right)} = {\frac{\quad{V_{0}\gamma^{2}}}{\quad{2\quad\pi\quad{\mathbb{i}}}}\quad{\oint{{q(z)}\quad{\mathbb{d}z}}}}} \\ {\quad{= {~~}{\frac{\quad{V_{0}\quad\gamma^{2}}}{\quad{2\quad\pi\quad{\mathbb{i}}}}\left\lbrack \quad{2\quad\pi\quad{\mathbb{i}}\quad\left( \quad{{{Residue}\left\lbrack {{q\quad(z)},z_{\quad 1}} \right\rbrack} +}\quad \right.} \right.}}} \\ \left. \left. {{{Residue}\quad\left\lbrack {{q\quad(z)},z_{2}} \right\rbrack} + 0 + 0} \right) \right\rbrack \\ {= {V_{0}{\gamma^{2}\left( {\left( \frac{\left\lbrack {z + \left( {{- 1} + a} \right)} \right\rbrack\left\lbrack {1 + {\left( {{- 1} + a} \right)z}} \right\rbrack}{\left( {z - z_{2}} \right)\left( {1 - {z_{1}z}} \right)\left( {1 - {z_{2}z}} \right)} \right)_{{|z} = z_{1}} +} \right.}}} \\ {\frac{\left\lbrack {z + \left( {{- 1} + a} \right)} \right\rbrack\left\lbrack {1 + {\left( {{- 1} + a} \right)z}} \right\rbrack}{\left( {z - z_{1}} \right)\left( {1 - {z_{1}z}} \right)\left( {1 - {z_{2}z}} \right)}\left. _{z = z_{2}} \right)} \\ {= {\frac{V_{0}\gamma^{2}}{1 - {z_{1}z_{2}}}\frac{\begin{matrix} {1 + {2\left( {z_{1} + z_{2}} \right)\left( {{- 1} + a} \right)} +} \\ {\left( {{- 1} + a} \right)^{2} + {z_{1}z_{2}} + {z_{1}{z_{2}\left( {{- 1} + a} \right)}^{2}}} \end{matrix}}{{- \left( {z_{1} + z_{2}} \right)^{2}} + {2z_{1}z_{2}} + {z_{1}^{2}z_{2}^{2}} + 1}}} \\ {= \frac{V_{0}\gamma^{2}}{\left( {N_{t}^{2}E_{c}\kappa\quad\gamma} \right)^{2}\left( {a - 1} \right)}} \\ {\frac{\left( {{2a} + {2N_{t}^{2}E_{c}\kappa\quad\gamma} - {3{aN}_{t}^{2}E_{c}\kappa\quad\gamma} + {a^{2}N_{t}^{2}E_{c}\kappa\quad\gamma}} \right)}{\left( {{2N_{t}^{2}E_{c}\kappa\quad\gamma} - {{aN}_{t}^{2}E_{c}\kappa\quad\gamma} - 4} \right)}.} \end{matrix} & ({A5}) \end{matrix}$ 

1. A method for pseudonoise code tracking, comprising the steps of: receiving a spread spectrum signal in a plurality of elements in a multiple element antenna array; despreading an output from each of the plurality of elements with an early pseudonoise code to form a plurality of early despread outputs; despreading an output from each of the plurality of elements with a late pseudonoise code to form a plurality of late despread outputs; accumulating the early despread outputs and the late despread outputs over a plurality of chip times; comparing an energy of the accumulated early despread outputs with an energy of the accumulated late despread outputs; controlling a chip time shifting amount based on the comparison.
 2. The method of claim 1, wherein the spread spectrum signal is received on a dedicated physical control channel.
 3. The method of claim 1, wherein the comparing step includes the step of calculating the difference of the square of the accumulated late despread outputs and the square of the accumulated early despread outputs.
 4. The method of claim 1, wherein the early pseudonoise code is advanced by approximately one half of a chip time with respect to an on-time pseudonoise code.
 5. The method of claim 1, wherein the late pseudonoise code is delayed by approximately one half of a chip time with respect to an on-time pseudonoise code.
 6. The method of claim 1, further comprising the step of weighting the output of each of the plurality of elements.
 7. The method of claim 1, wherein the plurality of elements comprises all of the elements in the multiple element antenna array.
 8. The method of claim 1, wherein each of the plurality of elements is separated from at least one neighboring element by approximately one half of a wavelength c the spread spectrum signal.
 9. An apparatus, comprising: a processor connectable to a multiple element antenna array, the processor being configured to perform the steps of: receiving a spread spectrum signal in a plurality of elements in a multiple element antenna array; despreading an output from each of the plurality of elements with an early pseudonoise code to form a plurality of early despread outputs; despreading an output from each of the plurality of elements with a late pseudonoise code to form a plurality of late despread outputs; accumulating the early despread outputs and the late despread outputs over a plurality of chip times; comparing an energy of the accumulated early despread outputs with an energy of the accumulated late despread outputs; and controlling a chip time shifting amount based on the comparison.
 10. The apparatus of claim 9, wherein the spread spectrum signal is received on a dedicated physical control channel.
 11. The apparatus of claim 9, wherein the comparing step includes the step of calculating the difference of the square of the accumulated late despread outputs and the square of the accumulated early despread outputs.
 12. The apparatus of claim 9, wherein the early pseudonoise code is advanced by approximately one half of a chip time with respect to an on-time pseudonoise code.
 13. The apparatus of claim 9, wherein the late pseudonoise code is delayed by approximately one half of a chip time with respect to an on-time pseudonoise code.
 14. The apparatus of claim 9, further comprising the step of weighting the output of each of the plurality of elements.
 15. The apparatus of claim 9, further comprising a multiple element antenna array.
 16. The apparatus of claim 9, wherein each of the plurality of elements is separated from at least one neighboring element by approximately one half of a wavelength of the spread spectrum signal.
 17. The apparatus of claim 9, further comprising a voltage controlled oscillator connected to the processor, wherein the chip time shifting amount is controlled by varying an input to the voltage controlled oscillator. 